
1 0 / 5>2 4 d 3 



SUBSTITUTE SPECIFICATION^'^ifi^ 



TITLE OF THE INVENTION 

METHOD AND ARRANGEMENT FOR DESIGNING A TECHNICAL SYSTEM 

CROSS REFERENCE TO RELATED APPLICATIONS 

[0001] This application is based on and hereby claims priority to PCT Application No. 
PCT/DE2003/002566 filed on July 30, 2003 and German Application No. 102 37 335.3 filed on 

S 

August 14, 2002, the contents of which are hereby incorporated by reference. 
BACKGROUND OF THE INVENTION 

[0002] The invention relates to a method and an arrangement for designing a technical system. 

[0003] In order to design a complex technical system it is often necessary to optimize the 
system with respect to a plurality of contradictory criteria. The criteria influence target functions 
of the system, such as, for example, manufacturing costs or efficiency. In addition, possible 
operating points of the system can be restricted by auxiliary conditions. This leads to the 
problem of determining a set of optimal operating points for the system, that is to say the set of 
possible operating points of the system with which it is not possible to optimize the operating 
points further simultaneously with regard to all criteria. From the set of optimal points, individual 
users can then select the most suitable operating points of the system for their applications 
while taking into account secret criteria or expert knowledge. 

[0004] A weighting method for optimizing technical systems with respect to a plurality of criteria 
is known from C. Hillermeier: "Nonlinear Multiobjective Optimization: A Generalized Homotopy 
Approach", Chapter 3.2, Birkhauser Veriag, 2001 ("the Hillermeier Chapter 3.2 reference"), 
wherein scaling parameters are employed to apply transformations to scalar-value optimization 
problems. This method has the disadvantage that it is numerically very involved, because very 
many scalar-value optimizations have to be performed. Furthermore, the selection and variation 
of the scaling parameters necessitates an interaction with a user and in this respect cannot be 
automated. 

[0005] A stochastic method for optimizing technical systems with respect to a plurality of 
criteria, wherein a stochastic differential equation is used to solve the optimization problem, is 
described in C. Hillermeier: "Nonlinear Multiobjective Optimization: A Generalized Homotopy 
Approach", Chapter 3,3, Birkhauser Veriag, 2001. This method has the disadvantage that it are 
very involved in numerical terms, because a multiplicity of quadratic optimization problems have 
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to be solved. A further disadvantage lies in the fact that with the method, the individual target 
functions are not weighted, as a result of which important information for selecting an optimal 
point is not available to the user. 

[0006] A homotopy method for optimizing technical systems with respect to multiple criteria, 
wherein in addition to weighting factors for the target functions, Lagrange multipliers are used in 
order to take auxiliary conditions into account, is known from C. Hillermeier: "A Generalized 
Homotopy Approach to Multiobjective Optimization", Journal of Optimization Theory and 
Application, Vol. 110/3, pp. 557-583, Plenum Press, New York, 2001 ("the Hillemermeier Vol. 
110/3 reference"). The disadvantage of this method lies in the fact that an interaction with the 
user is necessary and In this respect the method cannot be automated. 

SUMMARY OF THE INVENTION 

[0007] One possible object of the invention is therefore to create an automated and numerically 
efficient method for designing a technical system. . . 

[0008] The inventors propose a method for designing a technical system in which the technical 
system is modeled by a predetermined set of target functions which are dependent on 
parameters. In this modeling process each individual target function is weighted with a 
weighting factor. The method solves an equation system comprising the parameters and the 
weighting factors as variables in a variable space, with solutions of the equation system forming 
operating points of a solution space in the variable space. In the method the operating points 
are determined by a predictor-corrector method, according to which, starting from a first 
operating point, a predictor generated by a stochastic variable is determined in the variable 
space, and subsequently, in a corrector step, a second operating point is determined. The 
determined operating points are used here to design the technical system. The method can be 
used to design a new technical system or modify or, as the case may be, adapt an existing 
technical system. 

[0009] An advantage resides in the fact that the method is automated through the generation of 
the predictor by a stochastic variable and so there is no longer any need for intervention on the 
part of the user. The linking of the numerical predictor-corrector method with stochastic 
methods guarantees efficient use of the computer resources for calculating operating points of a 
technical system. 
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[0010] In an advantageous embodiment the predictor is predetermined by random numbers, so 
that in particular during the execution of the method a random number generator can be used 
and through this, the automation of the method is ensured in a simple manner. 

[0011] In a further particularly advantageous embodiment the random numbers are normally 
distributed. What this achieves is that the trajectory of operating points which forms in the 
solution space during the execution of the method ensures a uniform distribution in the entire 
solution space and so ensures good coverage of all possible operating points. As a result of the 
use of normally distributed random numbers, in particular a Brownian movement on the solution 
space can be modeled by the method. 

[0012] Preferably the operating points which are determined by the method are what are known 
as pareto-optimal points which cannot be optimized further in relation to all target functions. In 
the method, in particular the points with positive weighting factors in the solution space are 
selected as operating points. 

[0013] In a further advantageous embodiment the operating points must also satisfy one or 
more auxiliary conditions, with the or each auxiliary condition being represented by a further 
variable of the equation system in the variable space. In this case the auxiliary conditions can 
be equality auxiliary conditions and/or inequality auxiliary conditions. With inequality auxiliary 
conditions a slack variable is preferably introduced, by which the inequality auxiliary conditions 
can be transformed into equality auxiliary conditions. The use of slack variables will be 
explained in more detail in the detailed description of an exemplary embodiment. 

[0014] The solution space of the operating points is preferably a manifold, in particular a 
submanifold in the variable space. In the HillemermeierVol. 110/3 reference it is explained 
under what preconditions the solution space forms such a manifold. 

[0015] Since in particular at the start of the method a first operating point is present initially, in a 
special embodiment this first valid operating point is determined by a weighting method, the use 
of weighting methods already being known from the related art (see the Hillermeier Chapter 3.2 
reference). 

[0016] With the predictor-corrector method, which may be used in the invention, a tangential 
plane to the solution space is determined, preferably in the first operating point, and the 
predictor is then specified in the tangential plane. 
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[0017] In a development of the method, if a negative predictor with one or more negative 
weighting factors occurs, a new predictor Is determined by a reflection at a subplane of the 
solution space of the valid operating points. Through this, new regions of valid operating points 
can be determined, which operating points can be of particular relevance to the user in terms of 
secret supplementary criteria or his/her expert knowledge. 

[0018] In a preferred embodiment, in the reflection step a point of intersection of the trajectory 
that runs between the first operating point and the negative predictor with a subplane of the 
solution space Is determined. The tangential component of the vector spanned by the point of 
Intersection and the negative predictor to the relevant subplane of the solution space Is then 
determined, with those weighting factors which were negative for the negative predictor in the 
points of the subplane now being equal to zero. Next, the normal component, associated with 
the tangential component, of the vector spanned by the point of Intersection and the negative 
predictor is determined. Finally, the new predictor is determined by two times subtraction of the 
normal component from the negative predictor. 

[0019] A Newton method known from the related art, which method is easily convertible 
numerically, is preferably used for the corrector method. 

[0020] The operating points are preferably detemiined by iterations of the predictor-corrector 
method, with the second operating point of the preceding iteration step being used in a current 
iteration step as the first operating point of the predictor-corrector method. In this case the 
method is terminated by, for example, an abort condition. In an advantageous embodiment the 
abort condition is met when a predetermined number of operating points has been determined 
and/or a predetermined time limit has been reached. 

[0021] In addition to the above-described method for designing a technical system, the 
inventors propose an arrangement for designing a technical system by which the above- 
described method can be performed. In particular the method comprises a processor unit by 
which it is made possible for the predictor to be generated using a stochastic variable. 

[0022] The arrangement preferably comprises a random number generator for generating 
random numbers which represent the stochastic variable. 
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[0023] The inventors also propose a computer program product which has a storage medium 
on which is stored a computer program which is executable on a computer and executes the 
design method. 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0024] These and other objects and advantages of the present invention will become more 
apparent and more readily appreciated from the following description of the preferred 
embodiments, taken in conjunction with the accompanying drawings of which: 

Fig. 1 shows a flowchart of the method according to one embodiment of the invention for 
designing a technical system; 

Fig. 2 shows a diagram which illustrates the predictor-corrector method used in one 
embodiment of the invention; 

Fig. 3 shows a diagram which illustrates the reflection method used in an alternative 
embodiment of the invention, and 

Fig. 4 shows a processor unit for performing the method according to one embodiment 
of the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENT 
[0025] Reference will now be made in detail to the preferred embodiments of the present 
invention, examples of which are illustrated in the accompanying drawings, wherein like 
reference numerals refer to like elements throughout. 

[0026] A flowchart of a method for designing a technical system is shown in Fig. 1. 

[0027] First, in step 101, a description form of the technical system is selected. The description 
form comprises a predetermined number k of target functions f=(f1 ,...,fk), with each of the target 
functions being dependent on n predetermined parameters x^ to Xp of the technical system. 

The target functions are for example the investment costs f i and the efficiency f2 of the 

technical system. In this case the target functions are described by the following equation: 




Investment, cos ts{x) 



Efficiency(x) 



) 
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where x=(x,,..., x„). 

[0028] The parameters to Xn can be configuration parameters or operating parameters of the 
technical system. 

[0029] By the method the valid operating points which are used for the design of the technical 
system are determined by the optimization of the target functions with respect to the 
parameters, whereby not all target functions f^ to f^ can be optimized simultaneously since the 

optimization criteria are generally in competition with one another. 

[0030] The technical system is further limited in the valid operating points by a predetermined 
number m of auxiliary conditions h=(hi(x)...,hm (x)) which can be expressed by the following 
equation: 

h(x)=0 

where 0=(0,...,0) represents a zero vector. In this case what is involved is an equality 
auxiliary condition, with an inequality auxiliary condition also being able to be considered as an 
alternative. An inequality condition of this kind is, for example: 

h(x)<0 or h(x)>0 . 

[0031] In order to solve the optimization problems by inequality auxiliary conditions, a number 
m of slack variables s=(si,...,Sm) are introduced, by which the inequality auxiliary conditions can 
be transformed into the following equality auxiliary conditions: 

h(x)+s=0 or h(x)-s=0 

[0032] With the optimization method described in the present embodiment, the valid operating 
points are known as pareto-optimal points, which satisfy the following condition: 

min,,, = {X e y{"\h(x) = 0} 

[0033] It can be shown that the solutions of this optimization problem are the solutions of the 
following nonlinear equation systems: < 



h(x) 



= 0 
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[0034] In this case the auxiliary conditions are taken into account by the Lagrange multipliers 
A.=(Xi,...,A.m,) and the target functions fj are weighted with weighting factors aj, whereby care 
must be taken to ensure that the total of all weighting factors is normalized to one, 
i.e. Zty a, - 1 = 0 . In this case, in particular, the individual weighting factors aj can also be 

negative or equal to zero. The solutions of the optimization problem are therefore vectors 
(x,A,,a) in the (n+m+k)-dimensional variable space of the above equation system. 

[0035] As shown in the Hillemermeier Vol. 110/3 reference, under certain conditions the 
solutions of this equation system describe a (k-l)-dimensional submanifold M in the variable 
space. 

[0036] The below described numerical steps for determining valid operating points are 
essentially based on the homotopy method described in the Hillemermeier Vol. 110/3 reference, 
wherein a predictor-corrector method is used for determining pareto-optimal points. 

[0037] In step 102, proceeding from the description form 101 of the technical system, a first 
pareto-optimal point z is determined by a standard method such as, for example, the weighting 
method. 

[0038] In this first pareto-optimal point, in the next step 103, a (k-l)-dimensional tangential 
plane T^M to the manifold M of the valid operating points is determined in point z. Toward that 

end, a Jacobi matrix of the equation system F in point z is subjected to a QR factorizing. From 
this, an orthonormal basis { q1 ...qk-i } is then determined which spans the tangential plane. 

The individual numerical steps performed in this process are described in detail in the 
Hillemermeier Vol. 110/3 reference. 

[0039] In the next step 104, a predictor y is determined in this tangential plane, with the 
predictor - in contrast to the homotopy method described in the Hillemermeier Vol. 110/3 
reference - being generated by a normally distributed random number vector b of the dimension 
k-1 in the tangential plane. In this case the predictor y has the following form: 

y =z +(qi...qk-l)b 

[0040] Through the use of a random number vector such as this, a Brownian movement can be 
modeled on the submanifold M, with the Brownian movement being able to be represented 
approximately as follows: 
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dZ, ^€P{Z,)dB, 
where 

P(z) is a projection matrix onto the tangential plane T^M in the valid operating point z. 
e is a scaling factor, and 

B^jsSRl is a Brownian movement in the variable space. 

[0041] In order to model this movement, the k-1 -dimensional normal distribution N(Ok.-|, tAsl^.i) 
is selected for b, where the mean value 0|(.i is the (k-l)-dimensional zero vector and the 
variance is the (k-l)-dimensional identity matrix l^.-i multiplied by a step increment t^ of the 
Brownian movement and the scaling factor e. 

[0042] An alternative method of determining the predictor is first to determine a normally 
distributed random number vector in the (m+n+k)-dimensional variable space and then to 
project the vector into the (k-l)-dimensional tangential plane T2M. 

[0043] After this, in step 105, the predictor is projected with the aid of a corrector method, which 
is, for example, a numerical Newton method, onto the manifold of the pareto-optimal points. In 
this way a new valid operating point is determined on the manifold of the pareto-optimal points. 

[0044] The steps 103,104 and 105 are repeated iteratively, with the operating point determined 
in the preceding iteration step being used as the starting point for calculating a new valid 
operating point. 

[0045] In step 106 a check is made to determine whether an abort criterion has been met, in 
other words whether, for example, a predetermined number of iterations have been performed 
or a predetermined time limit has been reached. If this is not the case, a return is made to step 
103 and the next iteration is performed. This is continued until the abort criterion is met. 

[0046] Once the abort criterion has been met in step 106, in a next step 107 the set of 
determined pareto-optimal points is restricted to those points in which the weighting factors q 
are positive. 

[0047] From these pareto-optimal points, in a final step 108, the user selects an efficient 
operating point of the technical system appropriate to his/her requirements and the technical 
system is designed using this efficient operating point. 
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[0048] Fig. 2 shows a two-dimensional graphical representation of the predictor-corrector 
method, which may be used in thedesign method. 

[0049] In Fig/2, designates a pareto-optimal point on the submanifold M, with this point 
having been obtained in the i-th iteration step of the method. In order to determine a new 
pareto-optimal point, the tangential plane T^jM to the submanifold M is first determined in the 

point z>. The tangential plane is indicated by dashed lines in Fig. 2. In the next step 104, a 

predictor point y'"*"'' is then determined using normally distributed random numbers in the 
tangential plane T^jM. In the following corrector step 105, which can be, for example, a Newton 

method, the new pareto-optimal point z'"*''' is determined. The method is then continued, with 
the pareto-optimal point z^'^^ being used as the starting point for new predictor step. 

[0050] Fig. 3 relates to a variation of the method, wherein if predictors with negative weighting 
factors aj occur, a reflection is performed in order to determine a new predictor with positive a\. 

Fig. 3 shows this projection step being performed in a three-dimensional representation. 
[0051] Fig. 3 depicts a case in which, starting from a pareto-optimal point z, a predictor y^eg is 
determined which has a negative a\. This is illustrated graphically in that the section between 
the point z and the point y penetrates the tangential plane T^M in the point S. In this case the 
point S in turn lies on a subplane of the tangential plane T^M, for the points of which the 
coordinate aj has the value zero. In order to perform the reflection, the point of intersection S is 
determined first. This can be done using a projection operator which projects the aj component 

from a parameter representation of the straight line running through the points z and y. After the 
point S has been determined, the vector xneg running between S and y can now be determined. 
This vector is then dissected into the tangential component t to the subplane and into a normal 
component n. Thus, t = Xneg-n applies to the tangential component. The reflection step is then 

performed, with the new reflected vector Xneu having the same tangential component t as the 

old vector Xp^g and the normal component corresponding to the normal component n of the old 

vector Xpeg with the sign reversed. The new vector is therefore x^eu = t-n = (Xneg-n)-n = Xpeg 

- 2n. There thus results a new predictor ypeu which was reflected at the tangential plane T2M. 

yneu S Xpieu applies to the new predictor point ypeu- The above described reflection 
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method increases the numerical efficiency of the method since the generation of points with 
negative weighting factors aj is avoided and consequently the technical computing resources 
are used more efficiently. 

[0052] Fig. 4 shows a processor unit PRZE for performing the method. The processor unit 
PRZE comprises a processor CPU, a memory MEM and an input/output interface lOS which is 
used in a different way via an interface IPC: An output is made visible on a monitor MON via a 
graphical interface and/or output on a printer PRT. An input is made via a mouse MAS or a 
keyboard TAST. The processor unit PRZE also has a data bus BUS which provides the 
connection from a memory MEM, the processor CPU and the input/output interface lOS. 
Additional components such as, for example, additional memory, data storage (hard disk) or 
scanner can also be connected to the data bus BUS. 

[0053] The invention has been described in detail with particular reference to preferred 
embodiments thereof and examples, but it will be understood that variations and modifications 
can be effected within the spirit and scope of the invention covered by the claims which may 
include the phrase "at least one of A, B and C" or a similar phrase as an alternative expression 
that means one or more of A, B and C may be used, contrary to the holding in Superguide v. 
DIRECTV, 69 USPQ2d 1865 (Fed. Cir. 2004). 



